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We present an explicit observationally acceptable model 
for evolution from inflation to the present epoch under the 
assumption that the entropy and matter of the familiar uni- 
verse are from gravitational particle production at the end of 
inflation. This eliminates the problem of finding a satisfactory 
coupling of the inflaton and matter fields. Since the inflaton 
potential V(<j>) may be a monotonic function of the inflaton 
(f>, the inflaton energy could produce an observationally sig- 
nificant effective cosmological constant, as in quintessence. 

I. INTRODUCTION 

A satisfactory inflation model for the very early uni- 
verse has to account for the entropy in matter — the 
fields of present-day physics or their predecessors. The 
usual assumption is that the entropy comes from the de- 
cay of the inflaton field or fields <fi whose stress-energy 
tensor drives inflation. In this scenario the inflaton rolls 
toward a minimum of its potential where it oscillates, and 
the oscillations are damped by the production of quanta 
of fields coupled to (j>. The couplings have to be extremely 
weak, to protect the flat inflaton potential which is nec- 
essary for a successful inflation picture, but at the same 
time strong enough to allow efficient thermalization of 
the inflaton energy density. Models satisfying these cri- 
teria do not arise naturally in particle physics; they have 
to be constructed ad hoc. 

A less commonly discussed alternative is that the en- 
tropy in the matter fields comes from gravitational parti- 
cle production at the end of inflation.^] Since the inflaton 
does not interact with matter, and its energy density can 
roll monotonically toward zero, it is a candidate for a 
present-day effective cosmological constant. This picture 
merits broader discussion. To this end we present a spe- 
cific example that we believe satisfies the main observa- 
tional constraints. Our choice of the inflaton potential is 
presented in the next section. The results of a numerical 



In another terminology, matter fields end up in squeezed 
states after inflation [til]. 



integration of the evolution from inflation to the present 
are shown in Figure 1. 

Ford [^| first investigated gravitational particle pro- 
duction at the end of inflation, and pointed out that 
particles created this way could account for the entropy 
of the present universe. Inflation would be followed by 
a period of expansion dominated by the energy density 
p$ of the inflaton. Ford noted that if the kinetic part 
(f> 2 jl were dominant the energy density would decay with 
the expansion factor a(t) as p$ oc a -6 , so the expansion 
could become radiation-dominated before light element 
production. Spokoiny || added the idea that the inflaton 
could end up as the scalar field that Peebles and Ratra Q| 
proposed serves as an effective present-day cosmological 
constant. Spokoiny || presented examples of potentials 
V(<p) in which there is a transition from inflation to ki- 
netic energy-dominated expansion, and he considered the 
conditions under which the expansion later is dominated 
again by the inflaton. Since gravity waves are produced 
by inflation with the same energy density per component 
as the matter fields, the entropy at very high redshift 
must be distributed among enough matter field compo- 
nents that gravitational waves do not make a significant 
contribution to the expansion rate during light element 
production. We discuss this condition in §VI, along with 
the possible detection of the gravitational waves by fu- 
ture detectors. Other possible connections to observa- 
tions include the effect of the kinetic energy-dominated 
expansion on relict particle abundances, as discussed by 
Kamionkowski and Turner (|), and on baryogenesis, as 
considered by Joyce || and Joyce and Prokopec 0. 
There may also be a conceptual value to the postulate 
that the inflaton potential resides in a sector that is not 
coupled to matter fields, and so may assume a simple and 
possibly "natural" form. 

Terminology might be mentioned. Spokoiny || 
calls the kinetic energy-dominated period "deflationary;" 
Joyce H prefers "kination." Caldwell, Dave, and Stein- 
hardt H propose the name "quintessence" for a field that 
acts like Einstein's cosmological constant A. Since the 
model under discussion reduces the role of the inflaton 
to the essential operation of driving inflation, and adds 
the possibility that the inflaton ends up as quintessence, 
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we call the picture "quintessential inflation." 



III. KINETIC ENERGY-DOMINATED 
EXPANSION 



II. THE MODEL 

We consider inflation driven by a single scalar infla- 
ton field <j) that interacts only with gravity and itself by 
the potential V((f>). The matter fields are scalar, spinor 
and gauge fields of some grand unified theory, and are 
supposed to be in their vacuum states during inflation. 

The homogeneous part of the inflaton field satisfies 



dt 2 



a d<j) 



dV 



»)/ TO pi> 



a dt i 
where the expansion rate is 

H 2 = (a/a) 2 = ( P4> + p 
and the Planck mass is written as 

rripi = (8ttG/3)- 1/2 = 4.2 x 10 18 GeV, 
with H = 1 = c. The energy density in the inflaton is 

P0 = 4> 2 /2 + V(<f>), 

and p m (t) is the matter energy density after inflation. 
We adopt the inflaton potential 



(1) 

(2) 
(3) 
(4) 



V = A(0 4 + M 4 ) 
AM 8 



</> 4 + M 4 



for 4> < 0, 
for (b > 0, 



(5) 



At —<fi M this is a "chaotic" inflation potential |T(J; at 
4> 3> M it is a "quintessence" form of the type considered 
in Refs. @], §, §. We adopt 



A = 1 x 10 



-14 



(6) 



from the condition that present-day large-scale structure 
grows from quantum fluctuations frozen into (j) during 
inflation.^ The following sections establish bounds on the 
constant energy parameter M . The example in Figure 1 
assumes the present density parameter in matter is f2 m ~ 
0.3, with f2^ = 1 — 17 m ~ 0.7 in the inflaton. Then with 
the model for matter and radiation in §§V and VI, 



M = 8 x 10 GeV. 



(7) 



We assume the inflaton starts at <j> <C — m p \ and rolls 
toward zero. Inflation ends at (f> ~ —m p i when a sub- 
stantial part of the potential energy has turned into the 
kinetic energy of the inflaton. This is followed by a ki- 
netic energy-dominated peri od where V(6) and the mat- 
ter energy density terms in Eqs. (|]) to (^) are negligible. 



Here the solution to Eq. 
density is 



llf) is <j> oc a fi , and the energy 



2 /2 ~ Xm^ajaf, 



(8) 



where a x is the expansion factor at the end of inflation. 
This expression in the expansion rate Eq. (|^) yields 



1/3 



= 2 1/2 m p i \n(a/a x ) 



(9) 



The central point, as shown by Ford || and Spokoiny jj), 
is that the inflaton energy density decreases faster than 
that of radiation, so if the potential remains negligible 
the energy of the particles produced at the end of infla- 
tion eventu ally dominates. Our potential decreases only 
logarithmically with a, V ~ AAf 8 m~[ 4 [ln(a/a x )]~ 4 , so if 
the radiation remains subdominant V(<f>) approaches the 
kinetic energy density at scale factor a* given by 



V M J 



/mpi\ 
V M ) 



(10) 



The expansion then enters inflation from which it never 
recovers. Our model has a chance to worli^ only if the 
matter energy density p m (a x ) at the end of inflation is 
large enough to dominate the inflaton energy density at 
a = a r < a*. We now turn to the estimate of p m . 



IV. PARTICLE CREATION AND THE 
RADIATION ERA 

Particle creation at the end of inflation can be stud- 
ied by the standard methods of quantum field theory in 
curved spacetimc fH|| . For a massless scalar field \ de- 
scribed by the Lagrangian 



L = \(d,x?-\s,n x 2 , 



(ii) 



2 On scales of astronomical interest the only difference from 
"chaotic" inflation with V(cj>) = \<f) A is the increased expan- 
sion factor to the end of inflation. Since the adibatic density 
fluctuations from the frozen inflaton depend on the logarithm 
of the expansion factor jl^, the increased expansion is not 
significant. In conventional inflation the value of A from the 
COBE normalization jTlj] is three times the rounded value in 
Eq. (0). Again, the precise value of A is not important. 



3 One sees from Eqs. (g) and (£)]) that p$ after inflation re- 
mains dominated by the kinetic energy part until p m becomes 
important if the potential varies as rapidly as V oc e~^^ a 
with 4> c < m p i/(3v / 2). In our example the exponential vari- 
ation of V with cj> is replaced by the rapid functional change 
of V at (j) ~ because of the small value of M/m p i. 



2 



where 1Z is the scalar curvature, this has been done by 
Ford [^j in the limit of nearly conformal coupling, |£ — 
1 1 <C 1. He assumes that during inflation the spacetime 
is close to de Sitter, and that inflation is immediately 
followed by radiation-dominated expansion. His result 
for the energy density of the created particles is 



= RH*(a x /a)\ 



(12) 



where H x is the Hubble parameter, inflation ends at a 
a x , and the numerical factor isn 



i?~l(T 2 (l-6£) : 



(13) 



Ford's analysis [0 was extended to arbitrary values of £ 
[13 1 and to arbitrary power-law expansion after inflation 
[ 15 1 . In all cases R ~ 10 -2 . We assume minimal coupling, 
£ = 0, and adopt R = 0.01 per scalar field component. 

Free massless spinor and gauge fields are conformally- 
invariant, so their contribution to p m is suppressed rela- 
tive to scalar fields.^ Thus in a simple model of gravita- 
tional production of quanta with negligible rest mass the 
matter energy density is given by Eq. (U) with 



R ~ 10~ 2 N S 



(14) 



where N s is the number of scalar fields. 

We consider next the thermalization of p m . For A = 
10 -14 the Hubble parameter at the end of inflation is 



H x ~ A 1/2 m pl - 10 12 GeV, 



(15) 



at the earliest filled squares in Figure 1. The created par- 
ticles have typical energy e ~ H x (a x /a) and number den- 
sity n ~ Re 3 . As Spokoiny || noted, the particles inter- 
act by the exchange of gauge bosons and establish ther- 
mal equilibrium among the fermions and gauge bosons 
when the interaction rate na becomes comparable to the 
expansion rate H. Here, a ~ a 2 e~ 2 and a ~ 0.1 — 0.01 
is the gauge coupling constant. With H ~ H x (a x /a) 3 , 
thermalization occurs at 



ath/a x 



l R-^ 2 ~ (10 2 - 10 3 )iV s - 1 /2 ) 



(16) 



4 The "sudden" approximation, in which the de Sitter line el- 
ement is matched to the radiation-dominated one, is adequate 
for modes with wave numbers k <C H x and can be expected 
to give correct orders of magnitude for k ~ H x . Particle 
production is strongly suppressed for k ^> H x . The main 
contribution to p m is given by the modes with k ~ H x , and 
thus Eq. (^|) can be regarded only as an order-of-magnitude 
estimate. 

5 In the case of interacting fields, particle production can 
occur even if the field equations are conformally invariant. 
This is due to the conformal anomaly of the quantum theory 

& 



only a few orders of magnitude expansion from the end 
of inflation, at temperature 

Tth ~ pTM ~ R^aH x ~ 10 9 V 3 / 4 GeV. (17) 

During kinetic-energy dominated expansion the ratio 
of energy densities in matter and the inflaton is (Eqs. ^ 
and 111) 



Pml Pi, ~ \R{a/a x ) 



(18) 



Radiation-dominated expansion thus begins at expansion 
factor 



a r /a x — (XR) 



-1/2 



(19) 



at the second earliest filled squares in Figure 1. We 
require a r < a*, where a* is given by Eq.([l0|). With 
R ~ 0.01 — 1, this gives the condition 



M < (10 iJ - 10 i4 ) GeV, 



(20) 



well above the value for quintessence (Eq. |?]). The tem- 
perature at the start of radiation-dominated expansion 



is, from Eq. (14), 

T r - P U 4 (a r ) - Ai? 3/4 m pl - 10 3 iV 3/4 GeV. (21) 



V. EVOLUTION THROUGH THE RADIATION 
ERA TO THE PRESENT 



At a ~ a r the inflaton satisfies 

4> r ~ 2 1 / 2 m pl ln(Ai?)- 1 / 2 , 4> r ~ R 3 ' 2 \ 2 ml x . (22) 

In an extended period after a r the potential gradient term 
in the field equation ([!]) for <fi is negligible and the solu- 
tion for radiation-dominated expansion with the initial 
conditions (p)) is <j) = A - B{t r /t) 1/2 , with A « 4>r and 
B ~ m p i. That is, 4> is increasing quite slowly, at During 
matter-dominated expansion there is a stable increasing 
solution, 



= fcA 1 / 6 M 4 / 3 i 1 /3 j 



(23) 



where k = (72/5) in the radiation era and k = 3 1 / 3 
in the nonrelativistic matter era. The inflaton remains 
close to constant at <p ~ 4>r until this solution becomes 
comparable to <j> r , at 



t a ~ A-V2 m 3 M - 4 [ln(Ai?)- 1 / 2 ] 3 ~ 10 11 y, 



(24) 



for our parameters. Since this is much greater than the 
present expansion time is nearly constant from the end 
of kinetic-energy dominated expansion to the present, as 
in Figure 1. The present energy density in the inflaton 
thus is p$ ~ V(tfi r ), and the condition that this is com- 
parable to the total present energy density is 
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M ~ A- 1/8 m^ 4 if 1/4 [ln(Ai 1 , )- 1 ] 1/2 ~ 10 6 GeV, (25) 

close to the numerical result (Eq. f^). 

When the energy density p^, in the inflaton remains 
subdominant to a > a r the expansion factor from the 
end of inflation to the time when p^ reaches the nearly 
constant value V{4> r ) follows as in Eq. (|Io|), 



a»/a x ~ (m p i4> r ) ' /M 



2/3/ M 4/3 



(26) 



This is a*/a x ~ 10 18 for our parameters. 

It might be noted that the model is not sensitive to 
the power law index in the potential after inflation. The 
potential in Eq. (||) at cf> > may be generalized to 



V 



XM 4 



1 + {(j)/M) a ■ 



(27) 



Since V at a» < a < a Q is determined by the density pa- 
rameter Qj, = 1— fi TO = 0.7 and Hubble's constant, which 
fix the present value of p^,, the combination AM 4+Q / '</>" 
is independent of a. For a — 4 and A = f x 10~ 14 we 
require M/m p \ = 2 x 10~ 13 to agree with Q^. For the 
same value of A and a — 6 we need M/m p \ = 1.3 x 10 -10 , 
and, for a = 2, M/m p \ = 4 x 10~ 18 . For this range of 
values of the power law index a and the indicated scaling 
of M the evolution to the present is very similar to what 
is shown in Figure 1. 

At a > ao, when the inflaton energy dominates, the 
inflationary expansion can be described in the slow roll 
approximation, which neglects (f> in Eq. (|l]) and the ki- 
netic energy density term in Eq. (Q). This gives 



0«0 r (l + Ci/oO 1/4 , 



a (x exp 



C 



(1 + CH Q t) 



1/2 



(28) 



where C = 16m 2 , / '30 2 w 0.01 and we have used H a « 
V (4>r) I ' m pi for the present expansion rate. The expan- 
sion rate remains nearly constant and the growth of <fi is 
approximately linear in t until the time tb ~ l/CH a ~ 
10 12 y 

VI. GRAVITATIONAL WAVES 

Overproduction of gravitational waves is one of the 
dangers to be watched for in quintessential inflation. 
Gravitons are described by the same equation as a mini- 
mally coupled scalar field, so the energy density of gravi- 
tons at the end of inflation is twice that of a single scalar 
field because of the two graviton polarization states [jl6| . 
Quanta of the inflaton field (j> will also be produced, so the 
ratio of energy densities in the graviton-inflaton quanta 
and the matter density (Eqs. 
tion is 



12, p. 4ft at the end of infla- 
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(29) 



The annihilation of field quanta as the temperature falls 
below the mass conserves entropy (apart from the neu- 
trinos at z ~ 10 10 ), so the matter temperature varies as 
T m cx a^A/^ 1 / 3 , where Af(T) is the effective number of 
spin degrees of freedom in equilibrium at temperature T. 
Thus the mass fraction in graviton and inflaton quanta 
immediately before electron-positron annihilation is 

r„ = ^(a^a„)« A ^V /3 < 0.07. ( 30 ) 
Pm N s \NthJ 



At matter thermalization following the end of inflation, 
Ath ~ 10 2 — 10 3 . When light element production com- 
mences at z ~ 10 10 , N n — 10.75 (in radiation with 
M = 2, electron-positron pairs, and three families of neu- 
trinos). One extra neutrino family would increase the 
mass density by 16% . As indicated in Eq. (|30|), the 
standard model for the light elements requires that the 
gravitons add significantly less than this |l6| . 

In a minimal GUT, the only light particles with masses 
below H x are those of the standard model. Here N s = 4 
(the electroweak Higgs doublet), Nth ~ 10 2 , and r n ~ 
0.3, the equivalent of two new neutrino families. This 
is unacceptable within the standard model. A model 
with an intermediate symmetry breaking scale at energy 
;$ H x can give a substantially larger N s . Alternatively, 
in supersymmetric theories there is a complex scalar field 
for each chiral quark or lepton. For the minimal model 
this gives N s = 104 and r n ~ 0.01, well within the bound 
inEq. ©. 

The spectrum of gravitational waves differs from the 
usual inflation picture |l(| because of the long kinetic 
energy-dominated period. Gravitational waves are pro- 
duced with strain per logarithmic interval of frequency 
h g ~ H x /m p \ at Hubble length crossing during inflation. 
The amplitude is nearly constant^] until the proper wave- 
length re-crosses the Hubble length at expansion param- 
eter du, after which the amplitude decreases as l/a(t) to 
present value 



h g - (H x /mpi)au/a Q , H(a u ) ~ Lua /a u 



(31) 



at present expansion parameter a a and present proper 
frequency ui. Waves that pass the Hubble length at the 



6 More accurately, when the proper wavelength is much larger 
than the Hubble length, and a oc t 1 ^ 3 during kinetic energy- 
dominated expansion, the more rapidly growingsolution to 
the field equation h g + 3h g a/a = is h g — Int mJm. Here we 
ignore this slow evolution. A more careful analysis by Giovan- 
nini [J^| shows that it results in a logarithmic correction to 
the gravitational wave spectrum (|wj), h g oc lj~ 1//2 \n{uj x /io). 
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time of equality of mass densities in relativistic and non- 
relativistic matter have present frequency 



3 x 10 _16 fi m Hz 



(32) 



Waves that pass the Hubble length at the end of inflation 
and at the start of radiation-dominated expansion have 
present frequencies 



T Q R- 1/A ~ 10 11 Hz, 
XR 3/4 T ~ 1CT 3 Hz, 



(33) 



for our parameters. The part of the gravitational wave 
spectrum that starts oscillating during kinetic energy- 



(34) 



dominated expansion, at u> r < u> < u> Xl is |T 



Between u r and u eq the frequency dependence changes 
to h g tx uj^ 1 , the same as usual inflation models. 

The form of the gravitational wave spectrum at ui > to r 
can serve as an observational signature of quintessential 
inflation. The signal in our example is too small to be de- 
tected by interferometric graviational wave detectors un- 
der construction, including LIGO and VIRGO, but may 
be within reach of future detectors. 



VII. DISCUSSION 

The model meets several observational constraints. 
Primordial density fluctuations are essentially the same 
as in a chaotic inflation model pToj ] , and lead to the adi- 
abatic cold dark matter model for structure formation. 
This is viable but not yet shown to be unique. Perhaps 
topological defects that originate near the end of infla- 
tion, or the fields that produce the entropy in this model, 
also play a role in structure formation. The expansion 
becomes radiation-dominated at redshift z ~ 10 18 , well 
before light element production, and the energy density 
in gravitational waves may be tolerable. Baryogenesis 
remains an open issue, as in conventional inflation ||, 
j?J. The model is adjusted to present-day density pa- 
rameter OT = 0.3. There is considerable evidence for 
this low value |ft7j| . There is evidence for the zero space 
curvature required here and in usual inflation, from the 
redshift-magnitude relation for type la supernovae p8[ , 
p9| . This difficult measurement should be checked for 
consistency with other constraints on A. The reading of 
evidence from the anisotropy of the thermal cosmic back- 
ground radiation (the CBR) on open vs. cosmologically 
flat models still is mixed pQj, 

Our model fails to show the tracking of the inflaton 
and matter mass densities proposed by Zlatev, Wang and 
Steinhardt || . We have checked that in the range of val- 
ues 2 < a < 6 for the potential V a f a (Eq. W[) 



well after inflation the energy density in the infla- 
ton is quite close to constant over ten orders of magni- 
tude of expansion centered on the present epoch. During 
this time the inflaton pressure is quite close to — p$, so 
the field does not respond significantly to gravitational 
perturbations. Thus we suspect that, in contrast to the 
tracking case [p2[, it will be difficult to find classi- 
cal cosmological tests that distinguish our model from a 
Fricdmann-Lemaitre model with constant Lambda. 

There are arguably unnatural features of our model. 
First, the model parameters have to be tuned to ar- 
range inflaton domination at the present epoch. This 
problem is common to all cosmological scenarios with 
Q m < 1. Second, it is somewhat unnatural for a small 
mass M -C m p i to appear in the potential of the inflaton 
field </> interacting only with gravity. Both problems can 
be addressed in models with an extra inflaton field x i n 
which M is not a constant but a slowly varying function 
of x- Different parts of the universe would then ther- 
malize with different values of M, and a low value of M 
in our neighborhood could be selected by anthropic con- 
siderations. Regions with M ^> 10 6 GeV would become 
inflaton-dominated much earlier, and galaxy formation 
is such regions would be suppressed. This is similar to 
the anthropic selection of the cosmological constant, as 
discussed in p3[ . 

The quintessential inflation model simplifies the role of 
the inflaton by decoupling it from the matter. It remains 
to be seen whether this will aid the search for a believable 
physical basis for the inflaton. 
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FIG. 1. Evolution in quintessential inflation. The total 
mass density is plotted as a heavy solid line, the lighter line 
is the mass density in the inflaton, and the dashed line is the 
inflaton field value. The bottom curve is the density param- 
eter (the fractional contribution to the square of the Hubble 
parameter by the mass densities in relativistic and nonrel- 
ativistic matter). The end of inflation, marked by the ear- 
liest (left-hand) filled squares, is defined by the maximum 
of a(t)H(t) (the minimum value of the comoving Hubble 
length). Radiation-dominated expansion commences at the 
next squares. The next square indicates the epoch of equal 
mass densities in radiation and nonrelativistic matter. The 
last squares mark the present epoch at radiation temperature 
T ~ 3° K. The parameter M in the potential (Eqs. | to |) is 
chosen so the present value of the density parameter in matter 
is Qm = 0.3. The model for the matter assumes R = 1 (Eq 
|l4| ) and A/th — 1000 effective scalar fields at first thermaliza- 
tion. Annihilation of the extra fields at T — 1000 GeV causes 
the step in the mass density at redshift z ~ 10 16 (Eq. 
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